Accepted for publication in ApJL 

Preprint typeset using I^T^X style cmulatcapj v. 10/09/06 



O 
O 

(N 

C/3 



OH" 

I ' 

o 



o 
in 

o 
o 



X 



THE LUMINOUS AND CARBON-RICH SUPERNOVA 2006GZ: A DOUBLE DEGENERATE MERGER? 
M. Hicken 1 , P. M. Garnavich 2 , J. L. Prieto 3 , S. Blondin 1 , D. L. DePoy 3 , R. P. Kirshner 1 , J. Parrent 4 

Accepted for publication in ApJL 

ABSTRACT 

Spectra and light curves of SN 2006gz show the strongest signature of unburned carbon and one of 
the slowest fading light curves ever seen in a type la event (Ami5(i?) = 0.69 ± 0.04). The early-time 
Si II velocity is low, implying it was slowed by an envelope of unburned material. Our best estimate of 
the luminosity implies My — —19.74 and the production of ~ 1.2M Q of 56 Ni. This suggests a super- 
Chandrasekhar mass progenitor. A double degenerate merger is consistent with these observations. 
Subject headings: supernovae: general — supernovae: individual(SN 2006gz) 



1. INTRODUCTION 

The conventional picture for Type la supernovae (SN 
la) is that they are thermonuclear explosions of C-0 
white dwarf stars (WD). Two possible routes to explo- 
sion have been explored: single degenerate (SD) models 
in which the WD is nudged to explosion by accretion 
from a binary companion, and double degenerate (DD) 
models in which two WDs merge. 

SN la have been th e essential tool in establishing th e 
Hubble constant (e.g. JJha et al.lll999t iRiess et alj|2005j ). 
show i ng a sub-critical matte r density (|Garnavich et all 
Il998t iPerlmutter et alj Il998f) . and demonst rating that 
the cosmic expansion is acc elerating (e.g., IRiess et all 
11993 IPerlmutter et alJll999f >. Yet the precise nature of 
SN la progenitors remains uncertain. 

In the SD scenario, there is only a narrow range of 
physical conditions for accretion to increase the WD 
mass over time, (see lYoon fe Langerl I2003L and ref- 
erences therein). The DD pathway may lea d to neu- 
tron star formation (|Saio fe Nomotol [2004), but re- 
cent studies suggest explosion SN la is possible 

(|Yoon, Pod siadlow ski. fe Rosswog||2007f ). Both scenarios 
may contribute to the SN la population so it is important 
to consider their observable differences. 

The presence or absence of carbon and its velocity dis- 
tribution should be strong clues to the nature of the ex- 
plosion. For example, multi-dimensional calculations of 
pure deflagrations predict that unburned material should 
be m ixed into low- velocity layers (e.g., iGamezo et al.1 
2003). The fact that carbon is not commonly seen in 
deep layers suggests that most SN la are not a result of 
pure deflagration. 

Carbon is most often weak or absent in optical and in- 
frared spectra of SN la (Bra nch et al.1 .2003: M arion et al.l 
120061 iThomas et all 12006( 1. This suggests that burning 
to intermediate mass elements takes place in the outer 
layers of a C-0 WD for most SN la. The model that 
matches this feature best is a delayed detonation in a 
SD system since DD merger and pulsating delayed det- 
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onation (PDD) explosions are likley to leave significant 
unbu rned material in the outer layers (jKhokhlov et al.l 
I1993H . Conversely, the presence of carbon raises the pos- 
sibility of a DD merger. DD mergers could also have 
masses and luminosities that extend beyond the Chan- 
drasekhar limit. Two recent SN la may have been merger 
events: t he extremely l umino us SN 2003fg was inter- 
preted by iHowell et aD (|2006f l as arising from a super- 
Chandrasekhar (S-Ch) progenitor, and a compressed sil- 
icon layer in SN 2005hj provides indirect evidence fo r an 
unburned layer atop the silicon (Quimby et al. 20071). 

Here we present spectra and light curves 5 of the type la 
supernova 2006gz that begin two weeks before maximum 
light. They directly show the strongest signature of un- 
burned carbon ever reported and a low silicon velocity at 
very early times. The light curve is very broad and lu- 
minous, approaching the luminosity of SN 2003fg. These 
properties suggest that SN 2006gz may have arisen from 
the merger of two WDs. 

2. OBSERVATIONS 

SN 2006gz was discov ered bvlPuckett fe Pellonil (|2006f ) 
and independently by IWinslow fc Lil (|2006h on 2007 
Sept. 26.0 (UT). It was classified as a type la event by 
iPrieto et al.l (|2006al ) who noted an unusual double ab- 
sorption feature instead of the typical Si II 6355A line. 
The supernova was 12" west and 28" south of the center 
of the Scd spiral galaxy IC 1277, corresponding to a pro- 
jected distance of 14.4 kpc. It should be noted that the 
21-cm H I radio-based redshift of 0.0280 (|Havnes et al.l 
11997( 1 is incorrect. We report z=0.0237 ±0.0004 based 
on our optical spectrum of IC 1277, in close agreement 
with other optical redshifts reported in NED 6 . 

Seven spectra were obtained with the 2.4-m MDM tele- 
scope and The Ohio State University Boiler and Chivens 
CCD spectrograph beginning 2006 Sept. 28.1 (UT), 14 
days before maximum light. A 2" slit and a 150 line/mm 
grating were used. The resolution was ~15A. 

A series of 10 spectra was taken with the Fred L. Whip- 
ple Observatory (FL WO) 1.5-m Tillinghast telescope and 
FAST spectrograph (jFabricant et al.1 1998). also starting 
Sept. 28.1 (UT). A 3" slit and 300 line/mm grating were 
mounted. The resolution was ~7A. 

Beginning Sept. 29.1 (UT), 23 nights of UBVr'i' pho- 

5 Available at http: / /www.cfa.harvard.edu/supernova/SN2006gz 

6 http: / /nedwww. ipac.caltech.edu/ 
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tometry were aquired on the FLWO 1.2m telescope and 
Keplercam instrument. Eleven comparison stars were 
calibrated on 4 photometric nights in BV, 3 nights in 
r'i' , and 1 night in U , with a typical U-band uncertainty 
of 0.014 mag per star. Host galaxy subtraction images 
were taken in BVr'i' on 2007 May 27.4 (UT) and in U on 
June 27.4 (UT), after SN 2006gz had faded sufficiently. 

3. ANALYSIS 
3.1. Spectra 

The FLWO and MDM spectra taken at -14 days both 
show features typical of a type la supernova. As seen 
in Figure [U intermediate-m ass elements ( Si II, S II and 
Ca II) are clearly present. iPrieto et alj (|2006af ) noted 
that the most unusual aspect of the early spectra is the 
extra absorption seen redward of the Si II 6355A feature. 
Absorption features a t this wavelength have been identi- 
fied w ith C II 6580A (jBranch et al J 120031: iThomas et all 
l2006t ). Using the SYNOW code to generate synthetic 
spectra (jFisher et al.lll999h . C II and Na I reproduce this 
feature and the absorption at 5490A better than a com- 
bination of H I and He I. We adopt this interpretation. 
C II lines at 4745A and 7234A may also be seen be- 
fore -10 days. At two weeks before maximum light the 
C II feature has a rest-frame equivalent width (EW) of 
25A in SN 2006gz while in SN 1990N (jLeibundgut et all 
1991L Uefferv et alJll992t ) the feature is seen at 5A EW. 
Thomas et al.l ( 2006[ ) reports a strong C II detection 
in SN 2006D with 7A EW at -7 days. In the case of 
SN 2006gz, the C II absorption strength is comparable 
to the Si II EW of 37A at -14 days but quickly weak- 
ens and is undetectable past -10 days (Figure [2]). The 
Si II EW increases after the initial spectrum and reaches 
70A near maximum brightness. 

Unlike the luminous and slowly declining SN 199 IT, 
the spectra before maximum light are not dominated by 
Fe III lines. But the ratio of Si II 5972 to the 6355A line 
seen here does correspond to t he values seen in othe r 
very slowly declining SN la fe.g.. lGarnavich et al.ll2004| ). 
Except for the 5972A feature, many of the Si II lines 
in SN 2006gz are particularly deep and narrow, includ- 
ing 3858A and 4130A. There is also a strong Si III line 
(4560A) and Si III may also be the source of the absorp- 
tion observed at 5500A (jBranch et al.ll2003h . 

We measure the Si II velocity to be ^13100 kms -1 at 
-14 and -13 days, in contrast to ~15500 kms -1 in 
SN 1990N at -14 days (Figure HJ. But the SN 2006gz 
C II velocity is similar to the silic on velocity in SN 1990N 
and SN 1994D before -11 days (jPatat et al.lll996f ). The 
Si II velocity does not show the rapid decline that was ob- 
served in other SN la before -7 days. At -10 days it is at 
~12500 kms" 1 and declines by only 1000 kms" 1 to 11500 
kms -1 by maximum light. In four cases of SN la caught 
in early phases (1990N, 1994D, 2002bo, and 2003du) 
the Si II velocities at t < -10 days were 2000 to 4000 
kms"" 1 higher than that at maximum light (|Benetti et al.l 
120051 ). In SN 2006gz, the velocity of the Si II is unusu- 
ally low and constant at the earliest times, possibly as 
a result of overlying material for which the C II line is 
evidence. A tantalizing hint of an early plateau in ve- 
locity is suggested by the unchanging velocities in both 
Si II and C II between days -14 and -13. After -7 days, 
the Si II velocity slowly declines, with a slope similar to 
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Fig. 1.— SN 2006gz spectra at 5 epochs. SN 2003fg at +2 days 
and a SYNOW fit to the MDM spectrum at -14 days are shown. 
Dotted lines: C II AA4745, 6580, and 7234 blueshifted by 15500 
kins -1 . 
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Fig. 2.— Si II and C II velocities. Up until -12 days, Si II 
is suppressed in SN 2006gz while its C II velocity is comparable 
to that of Si II in SN 1994D and 1990N. SN 2006gz velocities 
are consistent within the error bars (not shown) for nights with 
measurements from both telescopes. 

SN 1990N and SN 1994D. 

I n the DP merge r mod els (and PDD models) studied 
by iKhokhlov et alJ (|1993f ) , the detonation of the dense 
C-0 occurs within a low-density C-0 envelope. The un- 
burned envelope compresses the outer, partially burned 
material, creating a density enhancement in the silicon 
layer which slows the photospher e from receding into 
lower- velocity layers. In SN 2005hi. [Quimbv et all (|2007f ) 
inferred the presence of overlying material from the ob- 
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Fig. 3.— SN 2006gz UBVr'i' light curves. B-band has a faster 
rise and slower decline than SN 1991T. Inset: B — V from +35 to 
+51 days, with E(B — V)o a i = 0.063 removed. The Lira relation, 
offset by E(B - V) host = 0.18 ± 0.05, is shown. 

served velocity plateau, but they did not have direct ev- 
idence for a carbon-rich layer. In SN 2006gz, carbon 
is definitely present, but the possible velocity plateau is 
very brief and very early. 

3.2. Light Curve, Reddening and Luminosity 

SN 2006gz has one of the broadest light curves (Fig- 
ure[3]) of any SN la ever measured: its 15-day _B-band de- 
cline from maximum light is Ami§(B) = 0.69±0.04, with 
t max {B) = 2454020.2±0.5 JD. The rise time of 16.6 ± 0.6 
days is very short for the slow decline and does not fit the 
positive correlation between rise and decline time in light 
curve fitters such as MLCS2k2 (iJha. Riess, fc Kirshneil 
12007ft and dml5 (jPrieto et al.ll2006bf ). IStrovinkl (|200"7f ) 
find evidence for a bimodal SN la rise time distribution 
based on a sample of 8 objects: 16.64 ± 0.21 days and 
18.81 ± 0.36 days. The decline rate is also outside the 
range used to construct these fitters. The i'-band peaks 
several days after B and does not appear to have the pro- 
nounced trough or second peak expected for a very lumi- 
nous SN la. For comparison, SN 2001ay has Amis(S)~ 
0.6 — 0.7 but its luminosity of My ~ —19.2 is significantly 
lower than expected for its decline ratc(Phillips et al. 
2002). Slow-declining and overluminous SN la are 
found almost exclusively in late-type galaxies while 
fast-declining and subluminou s SN la are usually 
found i n early-type hosts (e.g.. Ijha. Riess, fc Kirshnerl 
12007ft . iMannucci. Delia Valle. fc Panagial (|2006ft show 
that there may be 2 populations of SN la progenitors: 
a promptly-exploding component that dominates the su- 
pernova rate in star- forming galaxies, and a "tardy" one 
that gives rise to most of the SN la in older hosts. How 
the SD and DD pathways may relate to the correlation 
between luminosity and host type, and to the prompt 
and tardy progenitor components, needs further investi- 
gation. 

The luminosity and extinction of SN 2006gz in UBV 
are given in Table [TJ W e correct for Galactic redden- 
ing us ing the dust maps of lSchlegel. Finkbeiner. fc David 
(1998) and apply if-corrections. We use a Hubble-flow 
distance modulus of \x = 34.95 ± 0.04 (z cm b = 0.0234, 
H = 73kms" 1 Mpc- 1 , Q M = 0.3, Q A = 0.7, and a pe- 



culiar velocity uncertainty of 300 kms 1 ) to get abso- 
lute luminosities of M B = -19.17 ± 0.04, and M v = 
— 19.19 ± 0.04, before correction for host reddening. 

There are, however, 2 pieces of evidence for host red- 
dening. First, the Na I absorption at the SN position has 
an EW = 30 ± 0.15A, giving rise to E(B - V) host < 
0.15 ±0.08 (jTuratto. Benetti. fc Cappellarol l2002). Sec- 
ond, the B—V color shows normal, albeit slow, evolution. 
It increases by ~ 1.3 mag between maximum light and 
35 days after. We apply if-corrections and correct for 
time dilation and Galactic reddening in the B — V color 
curve. From +35 to +51 days, SN 2006gz has a colo r 
evolution consistent with the Lira relation (jLiral fl995ft . 
implying a host color excess of E(B — V)host = 0.18±0.05 
(Figure [3]). We adopt this while cautioning that it may 
bo incorrect if Lira's relation does not apply. Since the 
spectra are normal after -10 days, and the Lira-corrected 
B — V and U — B colors at m aximum are consi stent with 
those reported in Figure 9 of IJha et alJ IpOOfift . it is not 
unreasonable to assume that Lira's relation applies. 

The correct value of the reddening law for SN la is 
under active study (jConlev et al.ll2007ft . We use Ry — 
3.1 and Ry — 2.1 to give a wide range of plausible 
host extinctions and we calculate absolute magnitudes 
for SN 2006gz (see Table [TJ). Using R v = 3.1 gives 
M B = -19.91 ± 0.21, and M v = -19.74 ± 0.16, with 
B — V = —0.17 ± 0.05. This is comparable to the possi- 
bly S-Ch SN 2003fg, which had an absolute magnitude of 
M v = -19.85 ±0.06 (H = 73kms~ 1 Mpc- 1 ) and a color 
of_B — V^— 0.15, implying its host reddening was quite 
small. The case of R v = 2.1 gives M v = -19.56 ± 0.11 
for SN 2006gz. 

We also calculated a UVOIR light curve of SN 2006gz 
withLboi(tmax) = (2.18±0 39) 10 fe erg s" 1 for R v = 3.1. 
We follow the procedure of lStritzinger et aTl ([2006) to de- 
rive a 56 Ni mass (M Ni ) of 1.2O±O.28M . For R v = 2.1, 
Mm = 1.02 ± O.21M . The detonation of a 1.4M WD 
can produce - 0.9 - 1.0M o of 56 Ni (Khokhlov, A. 2007, 
private communication) . In figure [H we plot against 
Am 15 {B)for the 16 SN la in lStritzinger et al.l (l2006h and 
add S N 2006gz, SN 2003fg, and SN 2001ay ijPhillins et al.1 
12006ft . 

The DD merger models of iKhokhlov et all (|1993h con- 
sist of a 1.2 Mq WD detonation inside C-0 envelopes of 
different masses. They produce 0.63 M of 56 Ni. The 
envelope decelerates the layers of the exploded WD, in- 
creases their density, and thus raises the diffusion time. 
This produces broader light curves, as seen in SN 2006gz 
after maximum. Merger models with larger, S-Ch pri- 
mary WDs should be explored to see if they can ex- 
plain the properties of SN 2006gz. A different cause of 
the fast rise time and high luminosity could be an off- 
cen tered explosion, with the nickel- rich region towards 
us (|Hillebrandt. Sim, fc Ropkell2007ft . 

4. DISCUSSION AND CONCLUSION 

The three outstanding features of SN 2006gz are the 
presence of carbon, the low early-time silicon velocity, 
and the luminous, broad light curve. 

SN 2006gz shows the strongest evidence for unburned 
carbon of any SN la observed to date. The strong carbon 
feature has a higher expansion velocity than the silicon 
and fades quickly. The carbon is part of the photosphere 
at early times. The obvious inference is that the explo- 
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then it produced M Ni = 1.20 ± 0.28M Q (R v = 3.1) or 
Mm = 1-02 ± 0.21M Q (R v = 2.1). Combined with this, 
the presence of copious silicon and other non-iron-peak 
elements rules out a SD detonation and makes a S-Ch 
progenitor a good possibility. In both SD and DD cases, 
the centrifugal force of rotation can allow WDs to grow 
beyond 1.4M . However, the DD scenario is simpler in 
the sense that no helium or hydrogen shell burning is 
required to build up the primary WD. 

Taken separately, the three main features of SN 2006gz 
can be explained by various models. First, the unburned 
carbon is plausible to different degrees in the SD and DD 
scenarios. Second, the suppressed silicon velocity is pre- 
dicted by DD models. Third, the luminous and broad 
light curve with high implied 56 Ni mass is possible in SD 
and DD models. However, taking all of these features to- 
gether points to a consistent picture in which SN 2006gz 
is the result of the merger of two white dwarfs. 



Fig. 4.— 56 Ni mass v. Ami5(B)for both reddening cases o f 
SN 2006gz. Als o shown are 16 SN la from lStritzineer et aLld2006ft, 
and S N 2003fg HHowell et al.ll200B , and SN 2001av IPhillips et ah! 
120061) . SN 2006gz and SN 2003fg may come from S-Ch progenitors. 

sion did not burn all the carbon. This could be explained 
by the outer WD C-0 layer being accelerated and diluted 
during a deflagration phase. It is also possible that the 
envelope of the DD merger could be shocked and accel- 
erated within the first few seconds of explosion and give 
rise to the observed carbon. 

The unusually low and slowly declining silicon velocity 
at early times, combined with the deep and narrow silicon 
absorption features, have been predicted by DD merger 
models. 

The light curve is one of the broadest SN la ever seen 
with Am 15 (B) = 0.69 ±0.04. If it obeys the Lira relation 
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TABLE 1 

Extinction and Absolute Magnitude at Maximum Light with E(B - V)host = 0.18(05) 
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